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A systematic expansion of the Boltzmann equation for the diffusion of dilute tracers, in powers 
of the Knudsen number, carried out to next-to-leading order (NLO), gives the relativistic causal 
diffusion equation (Kelly's equation). Using dimensionless combinations of dynamical quantities, 
we show when the small NLO term plays a crucially important role. We proceed to show that a 
derivation of Kelly's equation from a microscopic theory of the correlation function of the number 
density of diffusers is possible. The correlator fulfils the Green-Kubo relation for the diffusion 
constant, as well as an f-sum which goes beyond a purely phenomenological causal theory. We argue 
that the construction generalizes to the full hydrodynamic equations. 

o : 

After Eckart [l[ and Landau cast Navier-Stokes hydrodynamics into manifestly Lorentz covariant form, it was 
O |. realized that there still remained an essential acausality in the theory 0, 0| • The problem was solved in principle by 
Israel and Stewart || who wrote down a phenomenological theory which has since been applied to heavy-ion collisions, 
the quintessential laboratory for ultra-relativistic hydrodynamics [f|. The hydrodynamic equations are obtained by 
1 making a derivative expansion of the transport equations. If the longest microscopic time scale in the problem is 
t-H \ a relaxation time tr, then the hydrodynamic equations are valid for times, r tr. When the Knudsen number, 
K = tr/t, goes to zero, the microscopic origin of the theory is effaced and what remains is a continuum theorv[To|. 
' (— i ' The usual equations of hydrodynamics are obtained in this way at leading order (LO) in K, i.e., at O(K ) Q. Here 
> . we discuss the version of the equations which appears at the next-to-leading order (NLO) in K. 

Two questions seem not to have been addressed fully in the literature, and are interesting enough that we consider 
them briefly here. The first question is the following: is the restoration of causality at NLO necessary, or is it an 
accident? The second question is: given a correlation function computed from the underlying microscopic theory, how 
does one show that the causal relativistic hydrodynamic theory is obtained at late times? We examine these questions 
here for the simplest of hydrodynamic equations — the diffusion equation, and its causal version: Kelly's equation. 
The answers are generalizable to all hydrodynamic equations, as we argue. 

Recall that the continuity equation for the density, n, of a dilute tracer inserted into a thcrmalized background 
fluid, 
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where J is the tracer current, expresses conservation of the number of tracer particles. Putting this together with the 
phcnomenologically established Fick's law, 

J(x,t) = -DVn(x,t), (2) 
where D is called the diffusion constant, gives rise to the diffusion equation, 

dn(x,t) 



dt 



DWn(x,t). (3) 



D is an example of a transport coefficient. A comparison of the strength of the two terms in the equation yields a 
dimensionless variable which we call the Wiener number, 

W =(£)/< DV2 "> = £ (4) 

where r is the time and L the distance over which the diffusion phenomenon is observed. For typical diffusive 
behaviour we expect W, to be finite, non-zero, and constant in time, which yields the familiar result L ~ V Dt. This 
violates relativistic causality for r < DW/c 2 , since L > ct in that case. This is shown in Figure Q] One suspects that 
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the diffusion equation is not valid at very early times, but there is nothing in the equation that prevents this. As is 
well-known, the local conservation equation (eq. [I} obeys relativistically causality. As a result, its violation in the eq. 
((3|) can be traced to Fick's law. 

Kelly Q examined the Boltzmann transport equation for the distribution function of a dilute gas of tracer particles 
in the background of a thermal gas in equilibrium. The transport theory is assumed to satisfy relativistic causality. 
The condition that the background gas has no flow selects out a special frame in which the gas is at rest. In the 
remainder of this note we work in this frame, giving up manifest Lorentz covariance for notational ease[llj]. The 
diluteness aproximation is used to retain only linear terms in the tracer density inside the collision integral. The 
collision term then yields up a natural microscopic scale of time through the collision frequency, v — 1/tr. 




FIG. 1: The parabola is a locus of constant W and contains the diffusive wake. The sets of straight lines are loci of constant M 
(the inner set for M — 1, and the outer for M — l/c s , i.e., the light cone) and denote the causal front. The lack of relativistic 
causality of the diffusion equation is seen in the fact that the parabola lies outside the light cone at small r. For any M, the 
front and wake can be distinguished at late times. 



It was shown that the expansion carried out to order K yielded the hyperbolic diffusion equation- 

d 2 n dn 



and, an expression for the diffusion constant, 



D = t r c s , 



(•5) 



(6) 



where c 2 is the RMS velocity of the tracer particles. Explicit solutions of the differential equation with this subsidiary 
condition show that relativistic causality is restored 

The relative importance of the causal (the second derivative in energy) and the frictional terms (the hrst derivative) 
is given by the dimensionless ratio 
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This is exactly as expected, since the causal term arises at NLO in K, and one obtains the diffusion equation by 
neglecting it. The relative importance of the causal term and the force is 
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Phenomena occurring at constant M 2 would include a causal diffusive front propagating with speed V = L/t. Thus 
M is the diffusive analogue of the Mach number. The three dimensionless constants are not independent, since, by 
construction, M 2 = KW. 
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The existence of a non-vanishing scale tr allows restoration of relativistic causality. On the light cone one has 
L = ct. The locus of constant W lies completely inside the light cone for a t* given by ^/WDt» < ct*. A little 
manipulation using the expression for D in eq. © brings this to the form = tr/t* < M 2 /W, where M c is the 
value of M when V — c and K r is the value of K at which causality is restored. One may note, with that for 
t < r* the frictional term in eq. ([5]) can be neglected, so that using eq. ([6]), the remaining terms become a wave- 
equation. Since c 2 < c 2 , relativistic causality is restored. More correctly, one notes that in the problematic region 
K > ~ 0(1), and hence the expansion breaks down. Thus, causality is restored by shielding the problematic 
region from the continuum equation. Investigation of r < r* must be carried out in transport theory, which is causal 
by construction. 

In Q it was noted that the solutions of eq. ([5]) consist of two distinct regions — the "wake", i.e., L oc y/r, where 
characteristic diffusive behaviour is seen, and the "front", i.e., L oc r, where causal behaviour is restored. This is 
precisely the result of our analysis here (see Figure [1]). The wake is characterised by the physics of fixed W, and 
the front by fixed M. W is, of course, inherent in the continuum theory. The appearance of a new, non- vanishing, 
dimensionless constant K implies the existence of a new constant KW, which leads to a shielding of the acausal 
region, as discussed above. This answers our first question: the restoration of causality is a necessary consequence of 
the NLO expansion in K. 

This analysis also tells us when to use Kelly's equation rather than the usual diffusion equation. Clearly, at 
late times, the differences between the two are minor in the region of the wake. The strength of the front decays 
exponentially in 1/K [4(, so this is also negligible at late times. At late times, therefore, the physics of the two 
equations is identical. Kelly's equation is to be used at early times, as the above analysis of causality shows. The two 
equations give different results at times when the front and the wake have not clearly separated. However, it is useful 
to be aware of the fact that at such times higher order effects, or the full Boltzmann equation, could be needed. 

We turn next to the second question and explore the relation between correlation functions and the hydrodynamic 
equation Ln = 0, where L is a differential operator in space-time. Although the computation is classic textbook 
material [|| , the connection with relativistic causal hydrodynamics has not been noted. The linear response function, 
i.e., the expectation value of the retarded density-density correlator, Xi is given by 

X (k,z) = x(0,0)[l+izG(k,z)}, (9) 

where z is a complex frequency in the upper-half plane, k a spatial momentum, and G(fc,z) is the Green's function, 
i.e., inverse of the operator L [8j. For the diffusion equation one has 

i , „ n iDk 2 X (0,0) 

Gfc,* = hence X (M) = , n L 2 ■ 10 

z + iDk 2 - z + iDk 2 - 

The spectral density is the limit as k — > of the imaginary part of x/k 2 as z approaches the real line from above. 
Hence, for the diffusion equation one finds that 

£>X(0,0) 

o{u) = . (11) 
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This is the basis for the Green-Kubo formula 



Dx(0,0)=£now(7(w). (12) 



Given any microscopic theory (for example, a quantum field theory specifying the interactions between the tracer 
particles and the constituents of the background fluid) which enables one to compute the spectral function, a(u>), the 
diffusion equation can be obtained if wcr(w) has a finite limit as uj vanishes. The Green-Kubo relation, along with the 
positivity of ua, provide a derivation of the diffusion equation from the microscopic theory. In a complete microscopic 
computation, the moments of the spectral density also exist. Since, the spectral density of the diffusion equation does 
not have well-defined moments, the diffusion equation does not care for these details. 

A straightforward treatment of Kelly's equation runs into subtleties. From eq. |(SJ) one sees that 
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(13) 



k 2 D 

Since the last factor diverges except for constant M 2 = truj 2 / Dk 2 , there seems to be no sensible spectral function 
otherwise. One could either proceed with caution along curves of constant M 2 , or use a regulator. 

The text-book regulator is the memory kernel method Q . One delays the response of a system to a disturbance by 



introducing the memory kernel — 

f* D 

J(x,t)= / dt'D(t-t'Wn(x,t'), D(t-t') = exp 

Jo T R 



f-t 

TR 



(14) 
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where the explicit expression for the memory kernel above constitutes the relaxation time ansatz. Putting this together 
with the continuity equation, one finds 

mi, ^ 1 (i, \ ^x(0,0)fc 2 J> x (0,0) 

G{k ' z) = z + l Dky(i- lZ r R y x(m) = z + 1 d(p- z V5)' ff(w) = a,(i + a,^)- (15) 

Note that this spectral function is perfectly regular as k — ► 0. Our earlier analysis of dimensionless combinations of 
variables is also reflected in this spectral function. To recover the diffusion equation, the spectral function in (eq. ITTj) 
must be accurate for arbitrarily small u>; so this is the limit K — > 0. The new term in the Green's function in eq. (|15|) 
is a correction to the diffusion equation to NLO in K = ujtr. Although the extra terms are parameterically small in 
K, they have important effects when M <C 1 is constant. These effects arise at oj oc k, which is precisely the region 
at which causal effects manifest themselves. 

Since uj<j(uj) is positive and the Green-Kubo formula is obeyed for M <C 1, this is an acceptable description 
of diffusion. The analysis is completed by noting that eq. (|14p implies that Fick's law must be replaced by the 
generalization 



d 

1 + T «di 



J(x,t) = -DVn(x,t), (16) 



which gives the NLO correction in K to the left hand side. Eliminating J between the continuity equation (eq. [T]) 
and the above relation then yields Kelly's equation (eq. [5]) , as noticed in Q . A phenomenological theory, such as eq. 
(JTHJ), leads to Kelly's equation. However, such a theory would have two material constants, D and r, whereas the 
microscopic theory has only one. 

While the Green-Kubo formula provides the correct reduction of the microscopic theory to the continuum, the 
spectral function in eq. (|15[) contains more. Its first moment exists, and (ui) = D/tr. Through Fourier transforms, 
the moments of the spectral function are related to time derivatives of the imaginary part of the linear response 
function, x"(x, t; x', i') = ([n(x, t),n(x' ,t')])/2, where n is the number operator in the microscopic theory. Each time 
derivative on n(x, t) becomes a commutator with the Hamiltonian. For the first moment, the double commutator is a 
generalization of the f-sum rule (also called the Thomas-Reiche-Kuhn sum rule) of atomic physics, and can be written 
as 

■wtr(w) = i \ rfx (HO, t), H(0, t)],n(x, t)}) . (17) 



x(o,o)J ^ 

A straightforward computation of the right hand side then leads to D/tr = (v 2 ), the RMS velocity, and reproduces 
eq. ^ for the diffusion constant. The f-sum rule knows about the microscopic constituents through the right hand 
side of eq. (|17|) . and hence can be used to constrain models. As an example, note that the fluctuation-dissipation 
theorem gives 2mT — D/tr, where T is the temperature and m the mass of the diffuser. For heavy quarks, such as 
the charm or bottom, one would expect that a computation of tr would immediately yield D through this relation. 
For diffusion of quantum numbers dominated by light quarks, such as the baryon number or electric charge, where 
the effective mass of the carriers is unknown, this relation potentially gives new information. Higher analogues of the 
sum rule exist at higher orders in K = wtr. Keeping successively higher moments, one can construct successively 
higher order terms in the analogue of Fick's law. This is a construction of effective dynamical theories at different 
time scales. Our second question is now answered. 

Finally we consider whether this analysis can be extended to full hydrodynamics. The major difference between 
the diffusion and Navier-Stokes equations is that the latter are non-linear. In the continuum limit, K — > 0, all 
hydrodynamic equations can be derived by combining three kinds of equations: continuity equations (such as eq. [T]) 
which relate a density and its current, constitutive equations (analogues of Fick's law) which relate the current to 
the density through a transport coefficient, and thermodynamic relations of various kinds, which serve to eliminate 
redundant variables. Non-linearities arise through the fact that the constitutive equations contain terms in powers 
of the velocities. However, the dissipative terms are constructed to be linear, and therefore the essential features 
discovered here in the diffusion equation generalize to hydrodynamics. 

Corrections to the continuum limit, in the form of NLO terms in K <C 1 can then be expected to yield results similar 
to those we have obtained for the diffusion equation Q. In particular, the memory kernel method is expected to give 
correction terms to the constitutive equations similar to eq. (|16p . As a result, the corrections to the hydrodynamic 
limit of the spectral function will be similar to those encountered in eq. (|15[) . All these arguments generalize, and 
the NLO corrections to the hydrodynamic equations are always expected to restore causality. In addition, the Green- 
Kubo formulae for the transport coefficients are not modified — a matter of some importance for the computation of 
transport coefficients from QCD using weak-coupling methods, lattice computations or the AdS/CFT correspondence. 
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The NLO corrections in K give rise to a connection between the relaxation time and transport coefficient through an 
f-sum rule. Since the f-sum rule constrains the microscopic degrees of freedom, they can be valuable in interpreting 
results obtained from lattice or AdS/CFT techniques. Such a systematic microscopic derivation of the hydrodynamic 
equations improves upon the Israel-Stewart formulation by relating various new material constants to the old ones 
through the f-sum rules. 

In conclusion, the systematic expansion of transport theory in the Knudsen number, K, which gives the usual 
hydrodynamic theory to LO in K, i.e., at 0(1), lacks relativisitic causality at early times [3, HI- This acausal region 
is necessarily shielded at next to leading order in K, i.e., at 0(K ) through the generation of a time scale, tr, and the 
breakdown of the equations for r ~ tr. Phenomenological theories can give the same causal hydrodynamic equations 
as the NLO expansion while introducing seemingly new material constants. The systematic expansion in K relates the 
new constants to the usual transport coefficients through f-sum rules. An order-by-order matching of any microscopic 
computation of a correlation function to the hydrodynamic equations is possible, and can be seen as the construction 
of successive effective long-time theories. 
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making the latter possible. 
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